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Abstract 

We apply the Darboux theory of integrability to polynomial ODE's of dimension 3. 
Using this theory and computer algebra, we study the existence of first integrals for the 
3-dimensional Lotka-Volterra systems with polynomial invariant algebraic solutions 
linear and quadratic and determine numerous cases of integrability. 

1 Introduction 

In this paper we are concerned with the 3- dimensional Lotka-Volterra [1,2] system (LV3), 
having first integrals formed with linear and quadratic polynomials. The LV3 is 

x = x(a\ + bux + b 12 y + b 13 z) = Pi(x, y, z), 

y = y(a 2 + b 2 ±x + b 22 y + b 23 z) = P 2 (x,y,z), (1.1) 
z = z(a 3 + b 3 ix + b 32 y + b 33 z) = P 3 (x, y, z), 

where the overdot means time derivative. This system typically model the time evolution 
of conflicting species in chemistry and biology [3], where the linear terms (also called 
Malthusian terms) denote growth (or decay) rates of each species independently of the 
others; the self-interactive terms, also called Verhulst terms (bu) represent the control 
on over-population of each of the respective species (such as cannibalism or depletion of 
resources), and the cross-interactive terms (bij,i ^ j) represent inter-species interactions 
(such as predator-prey). It has been extensively studied, starting with the pioneer works 
of Lotka [1] and Volterra [2], in the case where all the Verhulst terms vanish. Then in 
three dimensions the essential parameters reduce to only 6. A detailed study of this case 
can be found in Ref. 4 



There are many other natural phenomena modeled by (LI), such as the coupling of 
waves in laser physics [5], the evolution of electrons, ions and neutral species in plasma 
physics [6]. In hydrodynamics, they model the convective instability in the Benard problem 
[7]. Similarly, they appear in the interaction of gases in a background host medium [8]. 
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In the theory of partial differential equations they can be obtained as a discretized form 
of the Korteweg-de Vries equation [9]. They also play a role in such diverse topics of 
current interest as neural networks [10], biochemical reactions etc. The systems interest 
became crucial after the work of Brenig and Goriely [11,12] wherein they prove that a 
large class of ordinary differential equations implied in various fields of physics, biology, 
chemistry and economics, can be transformed into a Lotka-Volterra of greater order using 
a quasimonomial formalism. In the context of plasma physics, all the nonlinear terms 
represent binary interactions or model certain transport across the boundary of the system. 

In solving for the LV3, it is worthwhile to know, given a set of initial conditions, what 
its long-time asymptotic behaviour will be or whether stable periodic solutions exist. The 
existence of stable periodic orbits would be rather important for experimentalists wishing 
to obtain and maintain a stable oscillatory state. Since in general the solutions of the 
LV3 cannot be written in terms of elementary functions, the two questions of asymptotic 
behaviour and the existence of periodic orbits are rather hard to answer. Nor it is easy 
to explore the general solution using numerical schemes since one has to prescribe all 
the parameters of the LV3 in terms of real numbers. It will then be of most interest 
to possess constants of the motion, which, as such, contain the trajectories, and permits 
the localisation of the solutions or find the asymptotic state. Particularly, if a three- 
dimensional dynamical system admits a constant of motion, then the phase space is foliated 
into two-dimensional leaves and therefore certain types of irregular orbits cannot occur. In 
some cases it is even possible to find for a given three-dimensional dynamical system two 
functionally independent constants of the motion. Then the orbits are the intersections 
of two-dimensional flow invariant solutions and therefore nonchaotic. (Indeed, chaotic 
behaviour is associated with nonintegrability of the dynamical system.) In fact, obtaining 
a constant of motion corresponds to a partial integration and is interesting both from an 
analytical and numerical point of view. In this last point of view, obtaining a constant of 
motion is equivalent to reducing by one unit the dimension of the phase space. Since in 
nonlinear problems we are usually interested in a full exploration of initial conditions, any 
reduction of the dimension corresponds to a dramatic saving of numerical computation. 
Moreover, the existence of a constant of motion provides a welcome check of the numerical 
scheme with respect to its accuracy and stability. 

To introduce the adopted terminology, we must note that a constant of motion may 
be time-dependent. Usually this time-dependency appears in an exponential form [13,14]. 
Some authors [14] use the name of invariants, name which they apply regardless of the type 
of constant of motion. As here we are concerned only with time-independent constants 
of motion, we will keep for them the more usual name of first integrals. Much research 
effort has been devoted to the problem of finding time-independent and time-dependent 
constants of motion for hamiltonian and nonhamiltonian ordinary differential equations. 
The reason is that except for some simple cases, this problem is very hard and no general 
methods to solve it are known up to now. Nevertheless several approaches were developed 
in the last years, thanks to the use of computer algebra facilities. The most important 
among them being: specific ansatz for a invariant [13-16], singularity analysis [17,18], the 
Lie symmetry method [19,20], the linear compatibility analysis method [21,22], rescaling 
method [23,24] and the Darboux method [25-28]. Among the first mentioned methods 
one can cite the Carleman embedding procedure [13], the generalised Carleman [14] and 
the Hamiltonian methods [15,16]. Although they differ in the details of computation, all 
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these methods are based on an a priori hypothesis on the form of the invariant or the first 
integral H. 

In 1878 Darboux [25] showed how can be constructed the first integrals of planar poly- 
nomial ordinary differential equations possessing sufficient invariant algebraic curves. The 
Darboux method is based on the possibility of writing the invariant ( or at least an inte- 
grating factor) as the product of different algebraic functions /, raised at a given power \. 
It is on the form of these functions /, that it is introduced an Ansatz. The fi are deter- 
mined straightforwardly identifying to zero the coefficients of a polynomial expression in 
x,y,z . In that sense the Darboux method is not so different from some of the others cited 
above. Nevertheless the experiance shows that we have somehow divided the difficulties 
of the unavoidable identification. Typically the number of equations obtained is greater 
than the number of unknowns. In order to satisfy the whole system of algebraic equations 
we must then introduce some conditions on the parameters describing the system ( |1.1[ ). 
From a physical point of view the rule of the game is to find first integrals with as few 
conditions as possible on the parameters of the given system. 

Although the Darboux theory of integrability works for real or complex polynomial 
ordinary differential equations, we are concerned here with the existence of real first inte- 
grals of fll.lD when the parameters ai,bij of the system, the three dependent variables x, 
y, z, and the independent variable t (the time) are real. 

The paper is organized as follows. The main lines of the Darboux theory for three- 
dimensional polynomial differential systems are presented in Section 2. The first integrals 
constructed using the Darboux theory are build here exclusively with polynomial invariant 
algebraic solutions. In Section 3, the Proposition 1 exhibits the linear algebraic solutions 
of the LV3. The following sections contain the cases for which at least a first integral 
has been found. The cases where the LV3 has a single first integral are considered first, 
see Sections 4 and 5. In Section 4, the Theorems 2 and 3 concern the cases where all 
the fi are linear and in Section 5, the Theorems 4 and 5 are for the cases where f± is a 
quadratic polynomial, the first being when the conic passes through the origin and the 
second when not. In Section 6 we consider the cases of integrability, i.e. the cases where 
two first integrals coexist. Theorems 6-10 are when at most one of the first integrals 
contain quadratic algebraic solutions and Theorem 11 is for the cases where the two first 
integrals are formed with a quadratic algebraic solution. Finally in Section 7 we give our 
conclusions. 



2 Computational method 

The Darboux method is based on the existence of algebraic invariant solutions. Suppose 
that we can determine two polynomials fi(x,y,z) and Ki(x,y,z) such that 

+ + = m . (2.i) 

ox oy oz 

Then equation fi(x,y,z) = describes a surface which is formed by trajectories. It is an 
algebraic solution of the system. The polynomial K{ is called the cofactor of fi. 

These fi are going to be the "bricks" with which we will build the invariants. Suppose 
that we have obtained q functions /, . Let us consider the following function of the variables 
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x, y, z and t 



H = Y[f^(x,y,z)exp(st) . 



i=i 



Taking the total derivative 

dH dH dH dH dH 

—rr = -wr + -7-P1 + -7-P2 + -k-Ps 

at dt ox ay dz 



H 



ji v & x °y ® z 



(2.2) 



(2.3) 



Imposing that H is an invariant and making use of (^]l|), we obtain 



s + ^2 x i R i 



(2.4) 



»=i 



The equations in the Aj are now linear equations. How many fi do we need? The equation 
( |2.1| ) shows that if the system is of degree m, is at most of degree m — 1 independently 
of the degree of /j. So, the left hand side of (|2.4|) is a polynomial in x,y,z of degree 
at most m — 1 with a total of (m + 2)(m + l)m/6 terms. Consequently ( |2.4j) p roduces 
(m + 2)(m + l)m/6 equations where the unknowns are the Aj. We see from ( |2.4j ) that if 
we allow s to be different from zero the system of equations in the \ is inhomogeneous, 
and we need q = (m + 2)(m + l)m/6. If we want a first integral, then s = and the 
system becomes homogeneous and we need a priori either a new fi, or a new condition by 
imposing that the determinant of the system of equations in the Aj is zero. This will be 
the case here where we are concerned with first integrals exclusively. The required number 
of Aj in order to have first integrals is then q = (m + 2)(m + l)m/6 + 1. However, if it 
exists a solution with q = (m + 2)(m + l)m/6, then a first integral exists. In fact, all 
the first integrals presented here are with this value of q (Theorem 3 is for first integrals 
having a quadratic algebraic solutions which factorize to linear ones). We see that the 
possibility of solving for the A« in the system deduced identifying to zero the coefficients of 



the polynomial on x, y, z deduced from equation (|2.4! ) depends on the conditions we put 



on the coefficients of (1.1). We should not forget that some conditions come from equation 
( |2.4| ) and others from equation (|2.1|). 

We make a comment on the number of conditions on the system parameters in order 
to obtain a first integral and its variation with the degree of the invariant solutions tested. 
For the existence of invariant solutions of higher degree we require a larger number of 
conditions than for solutions of smaller degree. At least this is what happens for degrees 
one to three. For higher degrees one can get a saturation as this is what happens in 2d 
(see Ref 29). It turns out that the most interesting cases (i.e. the cases with a not too 
high number of conditions) will be obtained using straight lines and conies as invariant 
solutions for Lotka-Volterra and quadratic systems. 

In the case of system (|0|), m = 2 and in this paper we take q = 4. As the axes are alge- 
braic solutions, the problem of invariant search is reduced to the determination of only one 
algebraic solution i.e. f± = 0. Moreover f± is taken here as polynomial of degree at most 
two, namely fi = (l?-f) where 17* is a vector of dimension 10 and with components 1 or 
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and / a vector of components (/boo, hoox, fowV, fooiz, f20ox 2 , fnoxy, fwixz, f 020 y 2 , fonyz, 
foo2Z 2 )- 

The results presented here are given modulo the three-dimensional Lotka-Volterra sys- 
tem equivalences. This is because we can associate to a given LV3 ([O]) five equivalent 
three-dimensional Lotka-Volterra systems. The first two are obtained doing circular per- 
mutation of the variables x,y,z and of the parameters and bij, next three systems 
Vbij ^ are obtained doing the transformation 

(x, y, z, a 1 , a 2 , a 3 , bn, b 12 , b 13 , b 21 , b 22 , b 23 , b 31 , b 32 , b 33 ) -> 
(x, z, y, ax, a 3 , a 2 , bn, b\ 3 , b\ 2 , b 3 i, b 33 , b 32 , b 2 i, ^23, b 22 ), 

which keeps invariant the system, and the two others obtained doing circular permuta- 
tion of variables and parameters. We say that all these Lotka-Volterra systems are E 
equivalent. All the results of this paper are stated modulo these E equivalences. 

Table 1. Definition of the terminology used for the first integral conditions. 



na. 


definition 


na. 


definition 


ai 




lij 


bik — bjk — bkk 


Cij 


ai — aj 




bik + bjk — bkk 


c' ■ 


ai — 2a j 




2bik — bjk — 2bkk 


c" 


ai + aj 


i"' 

Hj 


2bik — bjk — bkk 


c"' 


ai + 2aj 


3ij 


bik + bjk — 2bkk 


Vijk 


ai + aj + a& 




bik + bjk - 46 fefe 


V 'ijk 


ai + aj — ak 


i" 


2bik + bjk - 4:b k k 


V ijk 


ai + aj — 2ak 


Hij 


bu(bij - bjj) - (b tj - 2bjj)(bji - ''„) 


u ijk 


ai + aj + 2a/c 


n'ij 


(bij ~ bjj)(b ki - b vl ) + (b^ + b k j - 2bjj)(bji - b vl ) 


Wijk 


2ai + aj + ak 


< 


b, ,(>>,, + b kj - 3b jj) - bki(kj - bjj) 


W 'ijk 


2a t + aj - a k 


< 


'■'•b.Ji,, + bjibkj - b u (bij + b kj ) - bjj(bji + b ki ) 


W ijk 


2a t - aj - a k 


'Hj 


(b i: j - bjj)(b k i - 2b,, \ - bjj(bji - bn) 


in'" 

w ijk 

bi 


ai + 2(aj + a k ) 


Oijk 


(b^ - bjj)(bj k - b k k)(bkt - bu) + (bji - bu)(b k j - bjj)(b ik - b k k) 


h 


Pij 


2bu(b k j - bjj) + bjj (bji - hi) 


dij 


bij ~ bjj 


V\j 


2bubjj - hj(b k i + hi) 


<j 


2bij — bjj 


V'lj 


b ki b k j - 2hj(2b k i - bji) 




bij — 2b jj 


v% 


bubkj — bki(bkj — bjj) 


e ij 


b ij + b jj 


Plj 


2hi(hj - bjj) + bjj(b k i - hi) 


e' 


2b tj + bjj 


Qij 


bijbubkk + bjibjj(hk — 2bkk) 


e" 


b^ + 2b,, 


q'ij 


hihk(hj — bjj) + bjjbjk(bji — bu) 


9ij 


hk - bjk 


q". 


bubjj(bkk + bjk) + 2b k i(bij - 2bjj) 


9ij 


hk — Zbjk 


a'.'! 


kibjj(bki - 2b vl ) - (b k j - 2bjj)(bji - 2bu) 


9% 


bik + bjk 


r ij 


Qibjj(bji - bu) + ajb il (b l j - bjj) 
ajbij — (ai + 2aj)bjj 


9% 


b ik + 2bj k 


r'- 

V 


hij 


bn bkj bjj bki 


r ". 


ajbij - (ai + aj)bjj 


Kj 


biibkj bjjbji 


r '.». 

V 


2a i b jl - (aj + aj)hi 




^j bij ^i bjj 


r iv 
ij 


aibjk - aj(b lk + 2bjk) 


K ij 


ajbij + aibjj 


Sij 


ajbij - (ai + aj + ak)bj 3 


hj 


aibjk — ajbik 


s' - 


2ajhj — (ai + 2aj + ak)bjj 




aibjk + ajbik 


s"- 


(ai + 2aj)bjj - aj(hj + b k j) 




a jbi] + "I, 1 ',, 


s"' 


(aibji — ajhi)(bkj — bjj)(bjk — bkk) + aihi(bjkbkj — bjjbkk) 
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3 Invariant algebraic solutions 

In this section we study the invariant algebraic solutions of the LV3 systems of degree 
at most 1. Thus, in Proposition 1, we present the invariant planes (invariant algebraic 
solutions of degree 1) together with their cofactors, and the conditions for their existence. 

Proposition 1. An LV3 has an invariant plane f = with cofactor K in the following 
cases, modulo the E equivalences, using the notation defined in Table 1 

(1) f = x = with K = b\\x + b\ 2 y + 613.2. 

(2) If b\2 = &13 = and a\bn 7^ 0, then f = a\ + b\\x = and K = b\\x. 

(3) If a\ = a 2 , 613 = &23 and d\ 2 d 2 \ 7^ then f = (6 2 i - hi)x - (6 i2 - 622)2/ = and 
K = ai + bux + b 22 y + h 3 z. 

(4) If a! = a 2 = a 3 , o 123 = and d 12 d 21 d 13 d 31 / 0, then f = (b 21 - b u )(b 31 - bu)x - 
O12 - b 22 )(b 31 - b u )y - (613 - 633X&21 - b u )z = and K = a 1 + b u x + b 22 y + b 33 z. 

(5) If r 12 = r 23 = r 31 = and a!a 2 a 3 bnb 22 b 33 ^ 0, then f = aia 2 a 3 + a 2 a 3 b n x + 
a\a 3 b 22 y + aia 2 &33-2 = and K = bux + b 22 y + 6332. 

(6) If b\ 3 = b 23 = r\ 2 = and aia 2 bnb 22 7^ 0, then f = a\a 2 + a 2 b\\x + a\b 22 y = and 
K = b\\x + b 22 y. 

(7) Ifa 3 = b 33 = r 12 = o U3 = and aia 2 b n b 22 ^ 0, then f = (6 31 -6n) (oia 2 + a 2 b n x - 
a\b 22 y) - a 2 b\\b\ 3 z = and K = b\\x + b 22 y. 

(8) If a 2 = b 22 = b 13 = b 23 = and ai&n&^cki 7^ 0, then f = (b 2l - &n)(ai + b n x) - 
bnbi 2 y = and K = b\\x. 

(9) If a 2 = a 3 = b 22 = 633 = 0123 = and 01611632^21 / 0, then f = 632(621 - &n)(ai + 
6nx) - 611612(6322/ - b 23 z) = and K = bux. 

(10) If a 2 = a 3 = 622 = 633 = 6 23 = 632 = and a^nc^i^i ¥= 0, then f = (b 21 - 
6n)(6 3 i - 6n)(ai + b n x) + 6n[6i 2 (6n - 6 3 i)y + 613(611 - 621)2;] = and K = b n x. 

(11) If a 2 = a 3 = 622 = 633 = 612 = 632 = and aibud 3 i 7^ ; then f = (6 3i - 6n)(ai + 
61121) — 611613Z = and K = b\\x. 

(12) If a 2 = 03 = 622 = 633 = 612 = d 21 = and 01611623632^31 0, then f = b 23 (b 31 - 
6n)(ai + 6hj;) + 6n6i 3 (6 3 2y - 6 23 ^) = and K = b n x. 

Proof. The proof is obtained finding both the linear / and K satisfying equation (2.1). 
So doing, we introduce in (|2.l| ) the general form for / and K, namely / = /boo + fioo x + 
/oioy + iboiz, K = kooo + kioox + kQioy + kooiz. The problem consits then in the evaluation 
of the coefficients of / and K by solving the algebraic system obtained identifying to zero 
the coefficients of the polynomial in x, y, z which results after this introduction of / and 
K in (2.1). The number of equations of this algebraic system is usely greater than the 
number of the unknown coefficients of / and K. In order to satisfy the entire system of 
algebraic equations, it is necessary then to introduce the conditions on the parameters of 
the differential system ( |1.1|) , which appear in the different statements of the theorem. ■ 
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4 First integrals formed with linear invariant algebraic so- 
lutions 

In this section we apply the Darboux method to deduce first integrals obtained with linear 
polynomial algebraic solutions. 

Theorem 2. Let fi = x, f 2 = y and = z. If a LV3 has a fourth algebraic solution 
/ 4 = of degree 1, then using the notation defined in Table 1, we can establish the following 
statements 

(1) // the LV3 satisfy to the conditions det (by) = and s = J2i=i a « a « = ®> where the cti 
are any non-trivial solution of (6y) T (aj) = 0, then the function \x\ OL1 \y\ a2 \z\ a ' i is a first 
integral. 

(2) If the conditions c± 2 = c 2 3 = 0123 = are satisfied and £^12^21^13^31 7^ 0, then the 
plane f 4 = (f>= (6 21 - 611) (631 - hi)x - (612 - 622X631 - hi)y - (613 - 633X621 - hi)z = 
is an algebraic solution of the LV3 and a first integral |/i| Ql |/2| a2 |/3| Q3 /4 exists, the cti 
being defined as solution of the system (b^) 7 \cti) = — diag(bij). 

(3) // the conditions c± 2 = g\2 = ^13 = are satisfied and di 2 d2i(|6ii| + I&21I + I&31I) 7^ 0, 
then the plane / 4 = (621 — b\\)x + (622 — b\ 2 )y = 0, is an algebraic solution of the LV3 
and a first integral \f 1 \ Xl \f 2 \ X2 \f3\ Xs \fa\ X4 exists with Ai = (b 21 - 611X613632 - 633622), A 2 = 
(612 - &22X613631 - 633611), A 3 = 613(612 - 6 22 )(6n -621), A 4 = 6i 3 [6 32 (6ii - 6 2 i) + 631(622 - 
612)] + 633(612621 - 611622). 

(4) If the conditions 613 = 623 = r i2 = are satisfied and a\a 2 b\\b 22 d\ 2 d 2 \ 7^ 0, then, 
the plane / 4 = 0102 + 026112; + a\b 22 y = 0, is an algebraic solution of the LV3 and a first 
integral |/i| Al |/ 2 | A2 / 3 A3 |/4| A4 exists with Ai = 6 22 (6 2 i -611), A 2 = 611(612-622), A 3 = 0,A 4 = 
611622 — 612621. 

(5) // the conditions r i2 = r 2 3 = r 3 i = 0103611622^23 + 0102622633^31 + 0203611633^12 = 
are satisfied and 010203611622633 7^ ; then the plane / 4 = 010203 + 02036112; + 0103622?/ + 
01026332 = is an algebraic solution of the LV3 and a first integral |/i| ai |/2| Q2 1/3|" 3 /4 
exists, the ai being defined as solution of the system (bij) T (a>i) = — diag(bij). 

(6) If the conditions 613 = 623 = di2 = d 2 \ = are satisfied and a\a 2 c\ 2 7^ 0, then, the 
plane / 4 = 0102 + 026112; + a\b 22 y = 0, is an algebraic solution of the LV3 and a first 
integral |/i| Al |/2| A2 /3 3 |/4| A4 exists with \\ = a 2 ,X 2 = — a±, A3 = 0, A 4 = a 4 — 02- 

(7) // the conditions a 2 = 613 = 622 = 623 = and 611612621^21 7^ 0, then the plane 
/ 4 = (621 — 6n) (ai + 6nx) — 6116122/ = 0, is an algebraic solution of the LV3 and a first 
integral / 4 Al 1/2(^/3 3 |/4| A4 exists with Ai = 0, A 2 = -611, A 3 = 0, A 4 = b 21 . 

(8) If the conditions bi 2 = 613 = 622 = 623 = are satisfied and aibnk 2 i 7^ 0, then 
the plane / 4 = a 4 + b\\x = 0, is an algebraic solution of the LV3 and a first integral 
|/i| Al |/2| A2 /3 3 |/4| A4 exists with Ai = a 2 6n, A 2 = -ai6n, A3 = 0, A 4 = 01621 - a 2 6n. 

(9) // the conditions c 32 = 612 = 613 = d 23 = d% 2 = ^23 = are satisfied and ai6n 7^ 0, 
then the plane / 4 = a 4 + b\\x = 0, is an algebraic solution of the LV3 and a first integral 
/i Al /2 A2 / 3 A; 7 4 4 exists with A 4 = 0, A 2 = -2, A 3 = 1, A 4 = 1. 

(10) // the conditions v' 231 = b± 2 = 613 = e23 = 632 = J23 = are satisfied and a\b\i 7^ 0, 
then the plane / 4 = a 4 + b\\x = 0, is an algebraic solution of the LV3 and a first integral 
fi x f 2 2 /3 3 / 4 4 exists with X 1 = 1, A 2 = -1, A 3 = -1, A 4 = 1. 
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(11) If the conditions c 2S = 612 = 613 = e23 = = J23 = are satisfied and a±bn 7^ 0, 
then the plane / 4 = a\ + b\\x = 0, is an algebraic solution of the LV3 and a first integral 
ft ft ft ft exists with A i = °> A 2 = -1, A 3 = -1, A 4 = 2. 

(12) If the conditions C23 = 612 = 613 = 623 = 633 = ^31 = O32 = are satisfied and 
^1^11^21 7^ 0, then the plane / 4 = a\ + b\\x = 0, is an algebraic solution of the LV3 and a 
first integral / Al |/2| A2 |/s| A3 |/4| A4 exists with X\ = 0, A2 = &n, A3 = — 611, A4 = 6n — 621. 

(13) If the conditions C31 = c 21 = 612 = 613 = d 2 \ = CZ31 = d' 32 = d 23 = ore satisfied and 
01611 7^ 0, i/ien i/ie pZane / 4 = a\ + b\\x = 0, is an algebraic solution of the LV3 and a 
first integral ft ft ft ft exists with A i = 0, A 2 = 1, A 3 = -2, A 4 = 1. 

(14) // the conditions c\2 = c' 3 \ = 613 = 623 = 633 = 631 = 632 = are satisfied and 
^12^21 7^ 0, then the plane / 4 = (621 — bn)x — (612 — b 22 )y = 0, is an algebraic solution of 
the LV3 and a first integral ft ft ft ft exists with A i = 0, A 2 = 0, A 3 = 1, A 4 = 1. 

(15) // the conditions c\2 = c 31 = b\\ = 613 = 622 = 623 = 631 = 633 = 513 = are satisfied 
and &12&21 7^ 0, then the plane / 4 = b 2 ±x — 612?/ = 0, is an algebraic solution of the LV3 
and a first integral ft ft ft ft exists with Ai = 1, A2 = 0, A3 = — 1, A 4 = 1. 



Proof. The cofactors of f\, f 2 , f% are K\ = a\ + b\\x + b\ 2 y + 613Z, K 2 = a 2 + 621X + 
+ 6232, K3 = 03 + bs\x + 6321/ + 6332; respectively. Statement (1) use only these three 
algebraic solutions. For the other statements, it is easy to check, that under suitable 
assumptions, / 4 = is an algebraic solution of the LV3, i.e. that / 4 verifies eq (2.1) 
where K4, the cofactor of / 4 takes the values = for statement (1) where / 4 = 1, 
K4 = a\ + b\\x + 6222/ + 633.2 for statement (2), K± = a\ + b\\x + b 22 y + b\%z for statements 
(3), (14) and (15), K4 = b\\x + 6222/ for statement (4), K4 = b\\x + 6222/ + 633Z for 
statement (5), K4 = b\\x for statements (7) — (13). We note that / 4 is the algebraic 
solution of Proposition 1(1) for statement (1), of Proposition 1(4) for statement (2), of 
Proposition 1(3) for statements (3), (14) and (15) of Proposition 1(6) for statements (4) 
and (6), of Proposition 1(5) for statement (5), of Proposition 1(8) for statement (7), of 
Proposition 1(2) for statements (8) — (13). Hence to each statement corresponds a first 
integral jx| Al |y| A2 \z\ Xi |/ 4 | A4 . ■ 



We must mention that statements (l)-(5) are, respectively, the invariants I, II, II', III' 
and III of Cairo and Feix [14] obtained using the Carleman method (with the additional 
condition s = for statements (1), (3) and (5). Note that statements (2) and (3) with 
a i = a 2 = a 3 = concerns the ABC system when 611 = 622 = 633 = (see for instance 
Ref. 22). With these additional conditions statement (2) concerns an integrable system 
(Theorem 6(22)). 

Theorem 3. The LV3 has a first integral \f x [ Al |/ 2 | A2 |/ 3 j A3 |/ 4 | A4 |/ 5 | A4 formed with the 
coordinate axes fi = x, f 2 = y and fj, = z and two other algebraic solutions / 4 = and 
fb = of degree 1 in the cases described by the following statements 

(1) // the conditions C23 = 611 = b 2 i = 623 = 631 = are satisfied and Ci3b 22 b 33 (a\b 22 b 33 + 
03(013(032 - b 22 ) - 612033)) / 0, then f 4 = a 3 + b 22 y and / 5 = (6 32 - b 22 )y + b 33 z and 
Ai = 03633(632 - 2622), A 2 = (632 - 26 2 2)(a 3 6i3 - 01633), A 3 = 01622633 - 03(612633 - 613622) 
and A 4 = -01622633 + 03(612633 + 6i 3 6 22 - 613632). 
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(2) If the conditions c 23 = 612 = &13 = 523 = s"^ = are satisfied and 01611^32^21 7^ 0, 
then f 4 = ai + b n x and f 5 = (b 32 - 622)2/ - (623 - 633)2 and Xi = 02(632 - b 22 )(b 23 ~~ ^33), 
A 2 = 01633(632 - 622), A3 = 01622(623 - 633), A 4 = 01(622633 - 623632)- 

(3) If the conditions 612 = 613 = ri2 = ^23 = r 3 ± = r' 2 " 3 = are satisfied and aia 2 a 3 c 23 bnb 22 b 33 
; then fi = a\ + b\\x and f§ = 010203 + 02036110; + 01036222/ + 01026332; and Ai = ; 
A 2 = 2a 3 /(a 2 - a 3 ), A 3 = -2a 2 /(a 2 - a 3 ) and A 4 = 1. 

(4) If the conditions C23 = 523 = h 23 = p[ 3 = ^12 = ^23 = are satisfied and 611612613622633 7^ 
0, then / 4 = 6 22 y + 6332 and / 5 = (ai + b u x)(b 2 i + 6n) + 6116222/ + 6116332; and Ai = 
-(621 + 6n)/6n, A 2 = A 3 = and A 4 = 1. 

(5) // the conditions C23 = 612 = 613 = 622 = d' 23 = g 23 = s 21 = are satisfied and 
01611623632 7^ 0, then / 4 = ai + b\\x and f 5 = 6322/ + 6232; and Ai = o 2 /ai, A 2 = -2, 
A3 = and A 4 = 1. 

(6) If the conditions a 2 = c 31 = b\ 2 = 613 = 622 = 631 = d2i = d 23 = are satisfied and 
0103611 7^ 0, then / 4 = a\ + 611X and /s = 01 + 6hj; + 26322/ + 6332: and Ai = A3 = 0, 
A2 = —2 and A 4 = 1. 

(7) // the conditions a 2 = a 3 = b\ 2 = 613 = 623 = 632 = d2i = 523 = are satisfied and 
01611622633 7^ 0, then / 4 = ai + b\\x and f§ = b 22 y — b 33 z and Ai = 0, A2 = A3 = —1 and 
A 4 = 1. 

(8) // the conditions a 2 = 03 = 612 = 613 = 633 = d' 32 = d 2 \ = <?23 = are satisfied and 
01611623632 7^ ; then / 4 = ai + b\\x and f$ = b 32 y + b 23 z and Xi = X 2 = 0, A3 = —2 and 
A 4 = 1. 

Proof. The first integral can be written as |/i| Al |/2| A2 |/3| A3 |/| A4 with / = / 4 / 5 . The 
cofactors of /1, f 2 , f 3 are Ki = a\ + b\\x + b\ 2 y + 6132;, K 2 = a 2 + 6 2 ix + b 22 y + 6232;, 
K 3 = 03 + 631X + 6322/ + 633Z respectively. Now it is easy to check that for statements 
(1), (2), (5)-(9) / 4 and /s are the algebraic solutions of Proposition 1 (2) and (3), for 
statement (3) / 4 and /s are the algebraic solutions of Proposition 1 (2) and (5), for 
statement (4) / 4 and /s are the algebraic solutions of Proposition 1 (3) and (5) and 
for statement (5) / 4 and /s are the algebraic solutions of Proposition 1 (2) and (7). 
Consequently / = is an algebraic solution of the LV3, i.e. / verifies eq ( |2.1| ) where 
K, the cofactor of / takes the value K = a 3 + 26nx + 2b 22 y + 26332; for statement (1), 
K = 03 + (621 + 6n )x + 6222/ + 6332; for statement (2), K = 2b\\x + b 22 y +2b 33 z for statement 
(3), K = ai(6 2 i + 6ii)/6ii + (62i + 6ii)x + 26222/ + 263 3 zfor statement (4), K = 2bux + 2b 22 y 
for statement (5), K = a 2 + b\\x + 2b 22 y for statement (6), K = 2b\\x + b 22 y + 633Z for 
statement (7) and K = 2b\\x + b 22 y for statement (8). For each statement, a solution 
of system Ai-fTi + X 2 K 2 + X 3 K 3 + A 4 iT = 0. Hence to each statement corresponds a first 
integral [x| Al \y\ X2 \z\ Xa |/ 4 | A4 \fa\ M . ■ 

The next two theorems concern LV3 systems having one first integral alone of type 
|y|A 2 |y 4 |A 4 y 4 q Uac i ra tic. These cases are characterized by n conditions 

among the parameters of the systems. The expression of all the corresponding first inte- 
grals are available in the e-mail address: "lcairo@labomath.univ-orleans.fr" . To prove each 
statement one must simply apply (|2.l| ) to obtain / 4 and its corresponding cofactors K 4 
and compute the Xi from the equation AiiTi + A2^2 + A3-K3 + A 4 ET 4 = 0, where K\,K 2 , K 3 
are the cofactors above mentioned for the coordinate axes, i.e. K\ = a\ + b\\x + b\ 2 y + 6132;, 
K 2 = a 2 + 6 2i x + 6222/ + 623^, K 3 = a 3 + b 31 x + b 32 y + 6332;. 
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5 First integrals with quadratic 

Theorem 4. Let fa = x, fa = y and fa = z. If a LV3 has a fourth algebraic solution, 
fa = of degree 2 passing through the origin, then a first integral |/i| Al |/2| A2 |/3| A3 |/4| A4 
can be obtained when the conditions given in Table 2 are fulfilled. 



Table 2: First integral conditions of Theorem 4. 
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Table 2: First integral conditions of Theorem 4 (continued). 
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( 3 ) di3d3i J23(6n 3 (2 632 613 + 6 22 633 - 3 633 632) - 6 22 6 3 i 3 di 3 - 6n 2 6 2 i (6 32 6 i3 + 2 6 22 6 33 - 3 6 33 6 32 ) + 
613 621 2 631 d 32 + 621631 2 di 3 rf 3 ' 2 + 6ii 2 6 3 i(3 6 2 2 633 + 2 6 33 632 - 4 6 2 2 613 - 632 613) - 611633621 2 c?32 + 
4 6n 622631 2 ^13 + 2 6n 621 631(2613 6 22 - 6 22 633 - 613 632)) ^ 

( 4 ) di3ji3(4 6n 6 23 d| 2 -6i2 2 6 2 3 d 3 i-6i2 6 23 6 3 i (3 632-4 6 22 )+6ii 6 i2 6 23 (3 6 3 2-4 622)+6i 2 6 i3 6 22 631- 
6n 613(2 6 32 2 - 632 6 22 - 2 6 22 2 ) - 6n 612 613 6 22 - 613 631^32(2 632 - 3 6 22 )) ^ 

( 5 ) &21&330Wl34'l4342521 ^ 
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( 7 ) aia 2 a 3 d'i 3 {b l2 d' 23 d' 2 ^ + 622(613633 - 2d| 3 )) ^ 

( 8 ) 6 3 34'i(«i-&234'i)^0 
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( n ) 623632(613632 + 3612633) ^ 

( 12 ) 611633^3 (&3 3 6i 2 + (46i 3 - 56 33 )6 2 2 ) ± 

( 13 ) 621633633(613 - 633^3) ^ 

Theorem 5. Let j\ = x, ji = y and fa = z. If a LV3 has a fourth algebraic solution, 
f\ = of degree 2 not passing through the origin, then a first integral |/i| Al |/2| A2 |/3| A3 |/4| A4 
can be obtained when the conditions given in Table 3 are fulfilled. 



Table 3: First integral conditions of Theorem 5. 
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Table 3: First integral conditions of Theorem 5 (continued). 
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Table 3: First integral conditions of Theorem 5 (continued). 
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6 Integrability 



We say that a 3-dimensional system like (IT) is integrable if it has two independent 



hrst integrals. Theorems 6-10 concern integrable LV3 having at least one first integral 
formed with linear algebraic solutions and Theorem 11 is for integrable cases with algebraic 
solutions being both quadratic. In these theorems we indicate the nature of the second 
first integral computed using the notation Pl(st) and Tk(st), respectively for Proposition 
1, statement st and Theorem k, statement st. and Q is for quadratic first integrals which 
do not belong to the statements of Theorems 4 and 5. 
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Theorem 6. An LV3 is integrable with two first integrals one of which being of the type 
of Theorem 2(1) in the cases of Table 4- 



Table 4: Integrability conditions of Theorem 6. 
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- &n = fei2 = &22 = &32 = <^23 = 523 


= 


633 


T2(6) 


(30) 


8 


as = 


c 12 ~ 


= &12 = 621 = b 3 i = b 3 2 = bs3 = g" 2 


= 


aia 2 


Q 


(31) 


8 


0-2 = 


r" - 
c 31 ~ 


'- bis = b22 = b 3 i = d'21 = d' 23 = g 3 \ 


= 


aia 3 


Q 


(32) 


8 


ai = 


6n = 


-- 621 = b 2 3 = b 3 i = dl[ 2 = di3 = r' 32 


= 


a2a 3 c 23 


Q 


(33) 


8 


ai = 


C23 = 


- fen = 621 = &31 = d\3 = ji3 = n 32 


= 


rf'2'3622633 


T4(2) 


(34) 


8 


a 2 = 


III 

c 13 = 


bi3 = 621 = b 22 = b 2 3 = 631 = 5i3 


= 


aia 3 


Q 


/OCX 

(35) 


o 

8 


a 2 = 


r" - 
c 13 ~ 


= &11 = 621 = b 22 = & 3 i = di 3 = g' 3 \ 


= 


7 

023512 


Pl(7) 


(36) 


8 


as = 


r" - 
c 21 ~ 


-- bi2 = 621 = b 3 i = & 32 = fe 33 = .921 


= 


(*) 


Pl(7) 


(37) 


8 


ai = 


as = 


bn = &13 = b 2 i = b 3 i = b 33 = d\2 - 


= 


622632 


T2(6) 


(38) 


8 


ai = 


as = 


bn = fei2 = &i3 = b 3 i = b 33 = d' 32 = 


= 


a 2 


Q 


(39) 


8 


a 2 = 


as = 


bl2 = &22 = b 2 3 = &32 = ^33 = 523 = 


= 


621 


Pl(8) 


(40) 


8 


a 2 = 


Oi = 


&22 = &21 = bl2 = fell = d' 23 = 512 = 


= 




Pl(10) 


(41) 


8 


ai = 


a 2 = 


hi = 612 = 621 = b 22 = b 3 i = du = 


= 


03512623613635 


Pl(8) 


(42) 


8 


a\ = 


a2 = 


03 = &33 = "31 = ^23 = n'23 = °123 


= 


612613623 


T2(2) 


(43) 


8 


ai = 


a 2 = 


a 3 = 612 = &22 = b 3 i = fe 3 2 = fe 33 = 





611 


T(2) 


(44) 


8 


ai = 


a2 = 


03 = bn = 621 = &22 = &23 = b 3 i = 





633 


PI (4) 


(45) 


8 


ai = 


a 2 = 


a 3 = &21 = &33 = d' 3 \ = n 23 = Oi 23 


= 


623 


Pl(4) 


(46) 


8 


a\ = 


a 2 = 


a 3 = &33 = d 32 = g'{ 2 = nsi = oi 23 


= 


611623 


Pl(4) 
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Table 4: Integrability conditions of Theorem 6 (continued). 



St. 


n 


conditions 


=U 














conditions 7^ 


f.i. type 


(47) 


o 
O 


ai = 


a 2 = 


03 = 


Oil = 


1~. 

013 = 


031 = 


033 = 


,7// n 

d 32 = 




L 

022 


T2(2) 


(48) 


o 

8 


ai = 


a 2 = 


O3 = 


611 = 


013 = 


l~ 

631 = 


632 = 


033 = 




7 

012 


Pl(4) 


(a n\ 

(49) 


o 
O 


ai = 


0-2 = 


O3 = 


O21 = 


1~. 

022 = 


023 = 


J 

"31 = 


J A 

"13 = 




331 


Q 


(50) 


Q 
O 


ai = 


a-2 = 


«3 = 


L 

"22 = 


033 = 


e3i = 


323 = 


r^f a 
n 23 = 




U h h s*?lt 

Ol2 t, 23 t, 32^3i 




(51) 


O 

O 


ai = 


a-2 = 


O3 = 


L 

Oil = 


i~. 

033 = 


in 

"32 — 


.923 = 


n 23 = 




T~.J-.T~. r,lt 

°13 23°225'l2 


T4(l) 


(52) 


O 

O 


ai = 


0.2 = 


O3 = 


013 = 


023 = 


«31 = 


ill 

«32 — 


n 12 = 




T~.1~.1~. 
Oll021&33 


Q 


(53) 


o 

8 


ai = 


a-2 = 


O3 = 


022 = 


033 = 


J// 

«31 = 


313 = 


/ A 

n 31 = 




777 

02lC»23"21 


Q 


(54) 


9 


Cl2 = 


j 

c 3 i ; 


7 

= 3 2 


T~ 

= 012 


1~ 

= 022 


= (121 


= Ol3 


= C(3i — 0S23 


= 


1 1 1 1 7 1 1 1 7 1 

|»l| + |on| + |o 33 | 


T2(5) 


(55) 





C12 = 


; C13 = 


- Oil 


u 

— 021 


031 


— O32 


— A" 

- "13 


jii „ 

— "12 — n 23 


A 

— U 


033022 




(56) 


9 


c 13 - 


- (-23 - 


= 6u 


= 621 


= 031 


= dl2 


= d 23 


= iT-32 = ^13 


= 


aia 2 a 3 


Q 


(57) 


9 


ai = 


r" - 
c 23 _ 


= &n = 


= &12 : 


= bi 3 = 


= b 2 3 = 


= b 33 = 


= e32 = .923 = 


= 




T2(l) 


(58) 


9 


a>3 = 




: 612 ; 


= 613 = 


= b 2 i = 


= b 23 = 


= b 3 i = 


= O32 = 633 = 


= 


|ai| + |&ii| + |6 22 | 


T2(4) 


(59) 


9 


as = 


c 2 'i = 


= 623 : 


= O33 ; 


= bi 3 = 


= d 32 


= ^31 


= d\2 = d 2 i 


= 


aia 2 


Q 


(60) 


9 


ai = 


r'" - 

c 23 _ 


= 611 = 


= 2 i = 


= b 23 ; 


= b 3 i = 


= 32 = 


= d\2 = d'{ 3 - 


= 


020-3 


Q 


(61) 


9 


a 2 = 


Ol = 


622 = 


= 2 1 = 


: b 32 = 


= &12 = 


: bi 3 = 


bn = d 23 = 





a 3 


Pl(8) 


(62) 


9 


ai = 


«3 = 


Oil = 


: 013 = 


= b 23 = 


; b 3 i = 


: b 33 = 


d\2 = d 3 2 = 





a 2 


Q 


(63) 


9 


ai = 


a 2 = 


a 3 = 


bn = 


022 = 


033 = 


023 = 


032 = 323 = 





I&21I + I&12I + 013 


T2(2) 



(*) |ai| + |6n| + I622I + |o 23 | ^ (**) |a 3 | + |6si| + |6 32 | + |6 33 | ^ 



Theorem 7. An is integrable with two first integrals one of which being of the type of 
Theorem 2(2) in the cases of Table 5, which exclude those already mentioned in Theorem 
6. 

Table 5. Integrability conditions of Theorem 7. 



St. 


n 


conditions 


=0 










conditions 7^ 


f.i. type 


(1) 


5 


Cl2 = C13 = 


= di2 


= J12 


= 0123 = 








Q 


(2) 


6 


C12 = C23 : 


= &12 


= 632 


= r 3 i = 0123 = 






011631633^23 


Q 


(3) 


6 


C12 = C13 = 


= O123 = d' 31 = d'21 = i" 32 = 






611633 


T4(49) 


(4) 


6 


C12 = C13 = 


= O33 


= n 'l2 


= n 23 = 0123 = 






d 3 i 


T2(4) 


(5) 


6 


&1 = 0-2 = 


O3 = 


.923 = 


h23 = °123 = 






622633di3d32 


T2(3) 


(6) 


7 


Cl2 = Cl3 = 


= &13 


= 6 23 


= o 3 i = o 32 = ri2 


= 




010203622612621 


Q 


(7) 


7 


ai = a 2 = 


03 = 


013 = 


023 = J23 = .931 = 







(*) 


T2(3) 


(8) 


7 


01 = a 2 = 


a 3 = 


&12 = 


013 = O23 = O123 = 


= 




633d2lj23 


T2(3) 


(9) 


7 


01 = a 2 = 


03 = 


022 = 


033 = 323 = O123 = 


= 




611623 


T2(3) 


(10) 


8 


Cl2 = Cl3 : 


= 011 


= b 22 


= o 3 i = d 2 ' 3 = 3 2 i 


= 331 


= 


621633 


T2(4) 


(11) 


8 


01 = a 2 = 


03 = 


&12 = 


di3 = d 3 i = d' 32 = 


p'1'3 = 





6iid 23 


T4(4) 


(12) 


8 


01 = a 2 = 


03 = 


&13 = 


2 2 = O23 = «32 = 


0123 = 





6 2 id 2 i 


T2(3) 


(13) 


8 


01 = a 2 = 


a 3 = 


Oil = 


032 = 033 = g 3 2 = 


0123 = 





622 


T2(3) 


(14) 


8 


01 = a 2 = 


a 3 = 


021 = 


di2 = d2 3 = e 3 i = 


.913 = 





633 


T2(3) 


(15) 


8 


01 = a 2 = 


03 = 


032 = 


033 =d' 21 = .923 = 


Ol23 = 


= 


6ll622di2 


T2(3) 


(16) 


8 


01 = a 2 = 


03 = 


022 = 


023 = d'21 = d 3 i = 


0123 = 


= 


633 


T2(3) 


(17) 


8 


01 = a 2 = 


03 = 


011 - 


°33 = d 3 ' 2 = 323 = 


0123 = 


= 


613623622 


T2(3) 


(18) 


8 


01 = a 2 = 


03 = 


022 = 


033 = d 31 = gi3 = 


0123 = 


= 


6nd 2 i 


T2(3) 


(19) 


9 


01 = a 2 = 


a 3 = 


Oil = 


b\2 = 033 = rf'3'2 = 


312 = 


323 = 




T2(3) 



(*) d 3 id32(6ii622 + 632d 3 'i) ^0 
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Theorem 8. An LV3 is integrable with two first integrals one of which being of the type 
of Theorem 2(3) (excluding those included in Theorem 7) in the cases of Table 6. 



Table 6: Integrability conditions of Theorem 8. 



St. 


n 


conditions =0 


conditions ^ 


f.i. type 


(1) 


5 


a x = a 2 = a 3 = 


.912 = 


i 1 " — 
'23 — 











^23jl3 


Q 


(2) 


6 


ai = c 23 = hi = 


= b 32 -- 


= ei2 = 


= .923 = 








&33^23^32^13jl2 


Q 


(3) 


6 


ai = c 2 3 = hi 


= &21 ; 


= &31 = 


= b 32 = 








(*) 


P1(2),T2(3) 


/ A \ 

(4) 


6 


oi = C23 = hi 


= 513 : 


= d' 12 


= 523 = 








111 J/f 

02a3«23"32"13"23 


Q 


(5) 


6 


ai = c 23 = hi 


= b 32 -- 


= 0123 


= 523 = 








^23^32^13 


Q 


(6) 


6 


ai = c 23 = hi = 


= &12 ; 


= 4 2 


= 523 = 








^13^23^12 


Q 


(7) 


6 


ai = c 23 = hi = 


= &13 ; 


= ^23 


= 523 = 








7 7 7//- 

h 2 di 2 d 32 j X3 


Q 


(8) 


6 


ai = c 23 = 6n = 


= 523 : 


= n 23 


= 0123 = 








b 22 b 33 d 23 d 32 d\ 2 d 32 


Q 


(9) 


6 


ai = c 23 = 6n = 


= h 2 -- 


= b 32 - 


= .923 = 








b 22 b 33 b\ 3 di 3 


T2(7) 


(10) 


6 


a 2 = C13 = fe 22 = 


= 523 : 


= 513 


= r 31 = 








b33 


Q 


(11) 


6 


ai = c 23 = 6n = 


= ^23 


= r 23 


= .923 = 








b 22 b 33 d\ 3 d 32 h' 23 


Q 


(12) 


6 


ai = c 23 = 6n = 


= &13 ; 


= b 23 - 


= .923 = 








b\ 2 d 32 


Pl(8) 


(13) 


6 


ai = c 23 = fen = 


= 523 : 


= 513 


= r 23 = 








a 2 a 3 d\ 3 d\ 2 


Q 


(14) 


6 


ai = a 2 = a 3 = 


512 = 


523 = 


531 = 








b 2 id 23 d 32 


Q 


(15) 


6 


Oi = a-2 = a 3 = 


&13 = 


&23 = 


523 = 










T2(3) 


(16) 


7 


ai = c 23 = fen = 


= b 32 -- 


= fel2 = 


= 512 = .923 : 


= 






b 33 b 2 \d 23 d 32 


Q 


(17) 


7 


a l = c 12 = c 13 ; 


= h' 31 


= P32 


= 532 = .931 


= 






b\ib 22 b 33 d\ 3 d\ 2 


T2(3) 


(18) 


7 


Oi = a 2 = a 3 = 


^32 = 


d'12 = 


"31 = 931 








( -T- -T- 1 


Q 


(19) 


7 


Oi = a 2 = a 3 = 


.923 = 


"23 = 


n 23 = Pl3 = 


= 






^12^13^22^23 


T4(l) 


(20) 


7 


Oi = a 2 = a 3 = 


&12 = 


&13 = 


n23 = 532 = 









111 Iff 

2 2033«32«32 


T4(l) 


(21) 


7 


ai = a 2 = a 3 = 


&12 = 


&32 = 


&33 = 523 = 









fe2lfel3 


T2(3) 


(22) 


7 


ai = a 2 = a 3 = 


4l = 


d' 31 = 


n 23 = 9 23 = 


: 






fe33^32d32*31 


Q 


(23) 


7 


ai = a 2 = a 3 = 


&13 = 


b 22 = 


d 23 = 523 = 











Q 


(24) 


7 


Oi = a 2 = a 3 = 


&33 = 


531 = 


.932 = 9 2 1 = 









7 7 J // 

0l3023Ol25l2 


Q 


(25) 


8 


C2 3 = C 31 = fen 


= b 21 


= b 3 i 


= b 32 = d' 12 


= di 3 


= 




010203^33^23 7^ 


P1(2),T2(3) 


(26) 


8 


Cl2 = C 23 = fell 


= r 2 i 


= r 31 


f f 

— r 23 = r 32 


= n'is 


= 




^22^33 


T4(37) 


(27) 


8 


Ol = c 23 = fen = 


= &13 ; 


= b 21 -- 


= b 3 i = d' 23 -- 


= 3i3 -- 


= 




020 3 


Q 


(28) 


8 


ai = a 2 — a 3 = 


&33 = 


e 3 i = 


e 3 2 = 513 = 


.923 = 







fe23^31&2lO'23 


T2(3) 


(29) 


8 


Oi = a 2 = a 3 = 


&13 = 


&22 = 


&23 = d' 3 \ = 


513 = 







b 2 i 


T2(3) 


(30) 


8 


ai = a 2 = a 3 = 


b 22 = 


a 21 — 


d' 13 = 512 = 


.923 = 







531 


T2(3) 


(31) 


8 


ai = a 2 = a 3 = 


hi = 


&22 = 


&33 = 523 = 


.913 = 







fli'2 


T2(3) 


(32) 


8 


ai = a 2 = a 3 = 


fell - 


&22 = 


&33 = 523 = 


n 23 = 







&23&21512 


Q 


(33) 


9 


c 12 = c 23 = b 12 


= b 22 


= fe 3 2 


— d 23 = e 3 i 


= .923 


= r 3 i 


= 


^33 


Q 


(34) 


9 


ai = c 23 = c' 3 \ 


= hi -- 


= hi -- 


= d'{ 2 = d 32 


= ^23 


= d 13 


= 


^21^33 


Q 


(35) 


9 


ai = a 2 = a 3 = 


&33 = 


d'i2 = 


d-21 = d 32 = 


.932 = 


9 f = 





hi 


Q 


(36) 


9 


ai = a 2 — a 3 = 


dl2 = 


dw = 


d 2 \ = d 3 i = 


: "32 = 


; Jl2 ~ 


= 


b 33 d\ 3 


Q 


(37) 


9 


a\ = a 2 = a 3 = 


b 22 = 


&23 = 


A" — A 1 - 
"21 — "13 — 


523 = 


5 3 'i = 





bi2 


Q 


(38) 


9 


ai = a 2 = a 3 = 


b 22 = 


&33 = 


d ' 3 \ = 523 = 


9& = 


9& = 





^12^23 


Q 



(*) fe22fe33O"23(022fel3 + fe33^12 - b 23 h 2 ) ^ (**) di 2 d 2 id 32 (h 2 b 21 - bnb 22 ) ^ 
(***) hidi3d'{ 3 i'{ 2 (2d 13 - b 33 ) ? (****) b 12 g' 3 \(b 12 b 21 + fenfe 32 ) ? 0. 



Integrability of the 3D Lotka-Volterra System 



527 



Theorem 9. An LV3 is integrable with two first integrals one of which being of the type 
of Theorem 2(4) in the cases of Table 7. 



Table 7: Integrability conditions of Theorem 9. 



St. 


n 


conditions 


Ij 
















conditions 7^ 


f.i. type 


(1) 


6 


bn 


= 612 = 


: b 2 i 


= fo31 


= fo32 = 


T23 


= 








a 2fo22fo33 




T2(7) 


(2) 


6 


bi3 


= 621 = 


- b 2 3 


= &31 


= ri2 = 


: 7-23 


= 








aia 2 a 3 




T2(4) 


(3) 


6 


bu 


= &21 = 


■- b 2 3 


= e 3 i 


= ^12 = 


: ^32 


= 








"2&11 




Q 


(4) 


6 


Cl2 


= c 32 = 


bi3 


= b 23 


= 0123 


= r w = 








aia 2 a 3 fond2i 




Q 


(5) 


6 


<Xi 


= C23 = 


foil 


= 923 


= 513 = 


^23 


= 








622^12 




Pl(3) 


(6) 


7 


C23 


= 6l2 = 


: fo31 


= b 32 


= e 2 i = 


: e 2 3 


= r 31 


= 






a 2 a 3 c'{ 2 




Q 


(7) 


7 


C31 


= ^12 = 


: fol3 


= b 3 i 


= &32 = 


e2i 


= r 2 3 


= 






aia 2 a 3 Ci 2 




Q 


(8) 


7 


c 23 


= ^23 = 


: fo31 


= b 32 


= fol2 = 


e2i 


= n 3 


= 






aia 2 a 3 Ci2 




Q 


(9) 


8 


r" 
c 12 


= 6ll = 


= &21 


= b 2 3 


= b 3 i = 


du 


= di3 


= r 32 


= 




aia 2 




Pl(2) 


(10) 


8 


r" 
C 31 


- foil - 


= &21 


= b 31 


= b 32 = 


du 


= du 


= »"23 


= 




«l"3fo33 




Pl(2) 


(11) 


8 


r" 

C 13 


-foil = 


: fol3 


= b 2 i 


= b 2 3 = 


&31 


= T 2 3 


= r'{ 2 


= 




aia 2 a 3 




Pl(2) 


(12) 


8 


C23 


= fol2 = 


: &22 


= b 3 i 


= b 32 = 


4-3 


= e 2 \ 


= s 13 


= 




fo23 




Pl(6) 


(13) 


8 


c 12 


= C13 = 


= fol2 


= b 2 i 


= b 2 3 = 


&32 


= e 3 i 


= ^13 


= 




aia 2 a 3 




T5(19) 


(14) 


8 




= C 32 = 


: fol2 


= fol3 


= b 32 = 


d'21 


— rl" 

— «23 


= d' 3 i 


= 




aia 2 a 3 




T5(14) 


(15) 


8 


C23 


— c 13 - 


: fol2 


= fol3 


= b 2 3 = 


d'21 


— rl" 

— «32 


= 932 


= 




aia 2 a 3 




Q 


(16) 


8 


C12 


— r'" — 

— c 23 - 


: fol2 


= fol3 


= b 2 \ = 


b23 


= b 32 


7/ 

= 4i 


= 




aia2d3 




T5(43) 


(17) 


8 


ai 


= foil = 


&21 


= &32 


= b 3 i = 


di3 


= d' 12 


= r 2 3 


= 




|"3| + |&33 




Pl(2) 


(18) 


8 


03 


= fol3 = 


&21 


= b 23 


= b 3 3 = 


d 3 2 


= e 3 i 


= r 12 


= 




a2 




Pl(2) 


(19) 


8 


03 


= fo33 = 


fol3 


= b 2 3 


= b 2 i = 


d 3 i 


— rl" 

— "32 


= r 12 


= 




aia 2 




Pl(2) 


(20) 


8 


ai 


= C23 = 


foil 


= b 2 i 


= b 3 i = 


di3 


= il3 


= r 32 ■■ 


= 




fo22 




Pl(3) 


(21) 


8 


02 


= fol2 = 


&21 


= b 2 2 


= b 3 i = 


&32 


= d 2 3 


= r 3 i 


= 




aia 3 




Pl(2) 


(22) 


8 


ai 


= foil = 


&21 


= b 3 i 


= b 32 = 


du 


= d'{ 2 


= S23 


= 




a2a3 




Pl(2) 


(23) 


9 


03 


= c 12 = 


fol2 


= fol3 


= b 2 3 = 


b3i 


= b 3 3 


= d' 21 = 


= d 32 = 





|ai| + |fon| + I&22I 


Q 


(24) 


9 


C23 


— c 31 = 


= fol2 


= fol3 


= b 2 3 = 


b32 


= b 3 3 


= d 31 


= r 12 = 





foll&22 




T2(5) 


(25) 


9 


r" 
c 12 


= C13 = 


= fon 


= b 2 i 


= b 2 3 = 


b3i 


= b 32 


= di2 


= d 13 = 


= 


aia 2 a 3 




Pl(2) 


(26) 


9 


Cl2 


= C23 = 


= fon 


= b 2 i 


= b 3 i = 


b32 


= d 12 


= di3 


Jff 

= "23 = 


= 


aia 2 a 3 fo22fo33 




Pl(2,3) 


(27) 


9 


Cl3 


= C23 = 


= fon 


= &21 


= b 2 3 = 


b3i 


= d 12 


= d'{ 3 


= ^32 = 


= 


(*) 




Pl(2,3) 


(28) 


9 


C23 


- C'21 = 


= fol2 


= b 2 2 


= b 3 i = 


b32 


= d 2 i 


= d 2 3 


— 3\2 ~ 


= 


|ai| + |fon| + 


633 


T4(4) 


(29) 


9 


Cl2 


= C13 = 


= fol2 


= b 2 i 


= b 2 2 = 


b 3 2 


= d 23 


= 4i 


= n 3 = 


= 


fon 




Q 


(30) 


9 


C21 


= C23 = 


= fol2 


= b 3 i 


= b 2 2 = 


b32 


= d 2 3 


= d 21 


= r 31 = 


- 


fol3 




Q 


(31) 


9 


03 


- r" — 

— c 21 — 


&23 


= b 2 i 


= b 3 3 = 


b\2 


= fol3 ; 


= d 32 ■■ 


- rl" 

- "31 — 





|"l| + |&ll| 




Pl(2) 


(32) 


9 


ai 


- r" - 

— c 23 — 


foil 


= fol3 


= b 2 i = 


b23 


= &31 : 


= b 32 = 
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Theorem 10. The LV3 has two first integrals with at least one being formed by one linear 
algebraic solution not included in Theorems 6-9 in the cases of Table 8. 

Table 8. Intcgrability conditions of Theorem 10. 



St. 


n 


conditions 


=u 














^ „ „ J __/ f\ 

cona. 7^ (j 


f.i. type 


(!) 


5 


ai = bu = 


7, 

612 


7. 

= o 32 


Jit 

= «13 










7, 7, 

03031033 


T2(7),Q 


(2) 


6 


L 7 

Oil = Oi2 = 


7 

= o 2 i 


7 

= 031 


= o 32 


/ 

— r 23 


= 






a203C 23 


T2(7),T5(l) 


(3) 


r> 




~ J 

a x = c 32 = 


7. 

011 


7 

= O13 


7 

= 23 


/ 

= .932 


= 






7 7 7 7// 

0l2022«12«32 


212(7) 


(4) 


7 


// 7 

C 23 = 023 = 


; 631 


= 32 


= O13 


= d 2 i 
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= 




aia2a 3 ci2 


T2(5),Q 
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; O13 


= 2 2 


= 631 


= b 32 


= (321 
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— «23 


= 


C12 


2P1(2) 


(6) 




C31 = O12 = 


; O13 


L 
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L 
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L 
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j 

= «3i 
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7 

t32 
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(7) 
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= O12 
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= 
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(14) 


8 
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= b 2 i 
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Theorem 11. The LV3 has two first integrals with quadratic algebraic solutions in the 
cases of Table 9. 

Table 9. Integrability conditions of Theorem 11. 
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n 


conditions 


=0 
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(1) 


7 


C23 = 012 = 


&13 


— «23 


— «32 


= 323 
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= 






(*) 


Q,Q 
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= d 31 
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611622 633 d23 
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7 


C23 = c'21 = 


012 
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d 32 d 31 
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7 
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= O31 
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= e 2 i 


= e 3 2 
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T5(30),Q 
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8 
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8 


r" — r" - 
c 12 — c 13 — 
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= ri2 
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= r 32 
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T4(2),T5(46) 


(8) 


8 
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= b 23 
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(9) 


8 
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(11) 
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(12) 


9 
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= r 32 


= 
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Q,Q 


(13) 


9 
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= n 23 
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T5(5), Q 


(14) 


9 


c 21 = c 31 = 
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(*) a 1 b 11 d 2 id 3 xd'i 3 d'i 2 + 0. 



We note, before concluding, that Theorem 4(6) and Theorem 8(1) with the additional 
conditions b\\ = 622 = °33 = and Theorem 8(31-32) are Theorem 6(5-6) and Theorem 
7(4) and 7(7) respectively of [4] in the case of ABC systems. Moreover, we note that to 



Integrability of the 3D Lotka-Volterra System 



529 



the above mentioned cases of Theorem 11, one must add the case obtained by Moulin- 
Ollagnier [30] for the ABC system which is formed with the first integral of Theorem 4(5) 
and an other one which is not of Darboux type (see Ref. 31 for more details). 

7 Conclusion 

The use of Darboux method for the LV3 has been rather efficient to find first integrals 
and to detect integrable cases. The main reason of this efficiency must be found in the 
fact that the problem of searching a first integral has been reduced here to the search 
of only one algebraic solution. (Theorem 3 is when a quadratic algebraic solution can be 
factorised in two linear ones). In total we have found 366 cases, among which 172 are cases 
of a single first integral and 194 when two first integrals coexist. As predicted, we see that 
for the existence of quadratic invariant solutions and the corresponding first integrals, we 
require a larger number of conditions than for linear algebraic solutions. In fact for having 
a single first integral formed by planes we can need only a maximum of 3 conditions in 
total (Theorem 2 (l)-(4)), whereas it is required 4 to 8 conditions when at least there is 
one second degree solution. The cases of integrability require at least 4 conditions when 
both first integrals are linear (in fact one case with 4 conditions and 5 with 5 conditions). 
Otherwise their number can grow now to 9 conditions. In general the number of conditions 
required is the number corresponding to the more constrained first integral increased by 
one. The number of essential parameters of the LV3 system being 9 (8 if we normalize the 
time), we omitted to consider the cases having more than 9 conditions. The present work 
can be completed with the cases where the planes x = 0, y = 0, z = and = are 
replaced by other planes and polynomials of greater degree, by exponential factors (see 
Cairo and Llibre [4]) and with time-dependent first integrals. 
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